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There has been a controversial recent claim that the standard result on the Higgs to two photon 
decay rate is incorrect, with the use of dimensional regularization fingered as the alleged culprit. 
Given the great importance of the ii" ^ 77 process as a possible Standard Model Higgs discovery 
channel at the LHC if the Higgs mass is light, it is critical to find a way to check the correctness of 
the results of dimensional regularization for this process. Here we report the results of a perturbative 
calculation of the // — > 77 decay amplitude using a spacetime lattice as a UV regulator, which is the 
only known gauge-invariant regulator for non-Abelian gauge theories other than dimensional regu- 
larization. We find that the decay amplitude calculated using lattice-regularized perturbation theory 
is consistent to very high statistical accuracy with the decay amplitude obtained using dimensional 
regularization. 



The decay of the Higgs boson into two photons, H — > 
77, is a crucial channel in the search for light Higgs 
bosons at the Large Hadron Collider (LHC) (see e.g. 
especially if the Higgs mass is close to the LEP bound 
of Mh > 114.4 GeV/c20. In the Standard Model (SM) 
there is no direct vertex for this decay, and the process 
proceeds through loops of intermediate gauge bosons and 
quarks. At one- loop level, the two most important con- 
tributions to the decay amplitude arc from the top quark 
loop and from boson loops. The evaluation of the 
latter contribution has recently been questioned. The 
M^^-loop contribution to the — >■ 77 decay amplitude 
takes the form 

where e and g are the U{1)em and SU(2)l couplings 
respectively, Mw is the VF-boson mass, kj^'^ are the four- 
momenta of the emitted photons, r = M^/4M^, and 
Mh is the mass of the Higgs boson. The function F{t) 
was first calculated in 0, HI many years ago, with the 
result 

Fir) ^l+r+^i2- r-i)[arcsin (3) 

for T < 1. These calculations used dimensional regu- 
larization, although [3l also derived some beautiful low- 
energy theorems on the properties of the amplitude which 
did not explicitly call on any particular regularization 
procedure. 

Recently this amplitude has been examined again by 
Gastmans, Wu, and Wu i. Rcf. i carried out the calcu- 
lation in unitary gauge directly in four dimensions with- 
out any use of dimensional regularization; they used a 
particular momentum-routing prescription to handle the 
divergences arising at intermediate stages in the calcula- 
tion of the amplitude. The result they obtained is 

^^M = ^ + ^(2-r-i)[arcsinV^]2, (3) 



which differs from the classic result by the absence of a 
term linear in t. This is a highly surprising result, since 
one expects that the calculation of a physical amplitude 
in quantum field theory must not depend on the regu- 
larization prescription. This difference was blamed on a 
problem with dimensional regularization. The claim was 
that dimensional regularization assigns a finite value to 
a contribution to the amplitude which is zero in d = 4 
when evaluated using the prescriptions of Q . If the strik- 
ing result of [sj were the correct form of the decay am- 
plitude, it would throw into question a vast number of 
calculations carried out with dimensional regularization. 
Furthermore, Eq. ^ gives a significantly smaller decay 
rate than Eq. ^ for Higgs masses where the diphoton 
decay is an important discovery channel. For example, 
for Mh = 115GeV/c^ the decay width would be reduced 
by « 50% 0. If correct, this would have significant im- 
plications for Higgs searches at the LHC. 

Given the long experience with the reliability of dimen- 
sional regularization, the claim of Q must naturally be 
viewed with skepticism, and a number of papers @ have 
appeared criticizing the results of Gastmans et al. from 
a variety of perspectives. However, given the extraor- 
dinary phenomenological importance of correctly deter- 
mining this amplitude in the Standard Model, it is highly 
desirable to add an extra check by explicitly computing 
using a reliable gauge- invariant regulator other than 
dimensional regularization. 

Lattice as a regulator — The list of reliable gauge- 
invariant regulators for non-Abclian gauge theories is 
very short. Naive momentum cutoffs break gauge in- 
variance, while Pauli-Villars-like regulators have seri- 
ous problems in non-Abelian theories (see e.g. 0]). 
In perturbative gauge theory calculations dimensional 
regularization is omnipresent since it manifestly pre- 
serves gauge invariance and is easy to use. The only 
other known gauge-invariant regulator is a spacetime lat- 
tice, which has the effect of imposing a momentum cutoff 
^ 1/a, where a is the lattice spacing, but does this in 
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a far more subtle way than a naive momentum cutoff, 
with manifest gauge invariance at any a. Lattice rcg- 
ularization has the great advantage that it can be used 
beyond perturbation theory, providing a nonperturbative 
definition of a theory. This has led to much progress in 
the nonperturbative understanding of gauge theories, e.g. 
using Monte Carlo methods to compute the full lattice- 
regularized Euclidean path integral of QCD. A lattice 
regulator can also be used in perturbative calculations, 
where it gives an alternative route to obtaining results 
that otherwise could only be obtained using dimensional 
regularization. 

Lattice regulators have two chief disadvantages for per- 
turbative calculations compared to dimensional regular- 
ization. The first disadvantage is that a lattice regula- 
tor breaks rotational and translation invariance, keeping 
only discrete subgroups of these symmetries. Rotational 
symmetry is crucial for enabling analytic evaluation of 
loop integrals in perturbative calculations; the loss of 
these symmetries on the lattice at finite a means that 
loop integrals must be evaluated numerically. This is 
certainly a practical inconvenience, but it is not a prob- 
lem of principle. The second disadvantage is of a deeper 
nature. It is not known how to implement lattice regula- 
tors for chiral gauge theories; for some recent reviews see 
0. Since the Standard Model is a chiral gauge theory 
this precludes the use of lattice regulators for computing 
most SM observables sensitive to the electroweak sector. 

Fortunately the difficulties with chiral fermion cou- 
plings do not play a role here. This is because only the 
W*-mediated contribution to the — >■ 77 decay am- 
plitude has been questioned, not the top-quark-mediatcd 
contribution. Using the lattice regulator we can sim- 
ply compute the _ff — > 77 decay amplitude with all the 
fermion fields turned off. 

Setup — We will carry out our calculations below the 
threshold for the H —J- W^W~ decay (r < 1) since 
the LHC has already excluded SM Higgs boson masses 
niH > 146 GeV[lO| at > 95% confidence. This as- 
sumption makes the Wick rotation to Euclidean space 
especially simple. (Above threshold, the calculation 
could also be done using lattice perturbation theory 
[m.) To evaluate the decay amplitude, we go to Eu- 
clidean space, choose fi = v ~ 1, and work in the Higgs 
boson rest frame, so that ki — (iA/^f/2, 0, 0, A////2), 
k2 = (iMjj/2,0,0, for on-shcll photons. The 
amplitude then takes the form 
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The lattice calculation proceeds by calculating F{t) 
numerically at a range of values of the lattice spacing a, 
taking the continuum limit Mwa, Mho, — > with r held 
fixed. With this end in view we calculate the dimension- 
less function F(t, aMw), where the dependence on aMw 
is due to lattice artifacts, and perform a two dimensional 



fit in r and aMyy to recover continuum limit, aM]y 0, 
giving F{t) = F{t, 0). For the fit function we take as our 
guide the functional form obtained from current analytic 
calculations and find it sufficient to consider 

F{t, aMw) ^Ci+C2T + 03(^2-^^ [arcsin ^/rf , (5) 

where ci, 02,03 are functions of aMw- 

For the gauge boson action that we use, a simple 
parametrization of the aMw dependence of ci{aMw) is 
found to be sufficient and C2 and C3 appear constant for 
the values of aAIw that we choose. 

Lattice Feynman rules — As in any perturbative cal- 
culation of an amplitude, we must fix a gauge to pro- 
ceed, and we choose to use unitary gauge. The lattice 
action for the gauge and Higgs sector of Standard Model 
in unitary gauge is described in Ref. [l2[. To carry out 
the calculation we use the HiPPy and HPSRC pack- 
ages developed for automated lattice perturbation the- 
We use the Symanzik-improved lattice ac- 



14 



ory I 

tion |15| for the vector-boson vertices that arise from the 
pure SU{2) gauge part of the action. The improvement is 
important since this greatly reduces lattice artifact con- 
tributions. These vertices are available in the automated 
packages and reduce to the continuum ones in the a — > 
limit. For the W boson propagator we add a mass term 
to the quadratic part of the gauge boson action and in- 
vert numerically to give the lattice version of the Proca 
propagator, as is appropriate to unitary gauge. 

The vector-boson mass term arises from the lattice La- 
grangian 



Lh = [av + aH{x + ae^)] [av + aH{x)] 

X Tr[eS 



(6) 



where v is the vacuum expectation value of the radial 
scalar field, H{x) is the Higgs field, Wi{x), i = 1,2,3 
are the SU{2)l gauge boson fields and B{x) is the U{1)y 
gauge field. The Ti are the anti-hermitian generators of 
SU (2) satisfying [Ti, Tj] = —CijkTk- As is usual for lattice 
actions Lh is dimensionless. The term quadratic in the 
W and B fields generates the and Z mass terms with 
Mw = Mz — gv/2 and where 



-cwZfj, + swAf^ . (7) 



Here sw = sin 0w , cw = cos 9w , where 9w is the Wein- 
berg angle, tan 9w = g' 1 9- This term also contains the 
usual HW^W~ vertex and, because the lattice spacing 
is non-zero, it additionally gives rise to extra interaction 
terms on the lattice which have no continuum counter- 
part. It should be emphasized that when inserted into 
a Feynman diagram these extra terms give contributions 
which are non-zero in the continuum limit; indeed, such 
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FIG. 1. Feynman rules derived from Lh (|6]). 




FIG. 2. Feynman diagrams (see text for description). 



contributions are a necessary consequence of maintain- 
ing gauge invariancc while regularizing the theory by 
restricting momenta to lie in the lattice Brillouin zone. 
These extra Feynman rules are derived by expanding the 
trace term in Lj^, using the Baker-Campbell-Hausdorff 
formula to the appropriate order and using ([7|) to identify 
the photon field A^. The HIV^W' and additional Feyn- 
man rules are displayed graphically in Fig. [1] The rule in 
(d) arises because of the presence of H{x + ae^) in Lh, 
and, because the Higgs boson is at rest, the Lorentz in- 
dex on the M^-bosons attached to this vertex is restricted 
to /.J = as shown in the figure. 

Summary of diagrams — In our unitary-gauge lattice- 
regulated calculation there are seven Feynman diagrams 
for the decay process, shown in Fig. [2] In the triangle (a) 
and turnip (b) diagrams, the vector boson vertices arise 
from the pure SU{2) gauge vector boson action and are 
similar to ones that arise in QCD. The ankh (c) and lat- 
tice turnip (d) diagrams arise from the additional lattice 
Feynman rules. Diagrams (e), (f) and (g) are modifica- 
tions of (a), (b) and (d) to include the lattice forward 
derivative coupling due to the (Vl^'^ H)W'^W~ term 
in the action. Note that there is no such modification 



of the ankh diagram as it would require the photons to 
have Lorentz index /x = 0, which they do not. 

Cancellation of divergences — The result of combin- 
ing these diagrams must be finite since there is no tree- 
level process for H 77, so no counter-term can be 
constructed to cancel any remaining divergences through 
the usual rcnormalization prescription. Therefore diver- 
gences {aMw)~"' present in individual diagrams must 
cancel. To investigate this cancellation in more detail, 
we fix r and study the dependence on aMw. By verify- 
ing the cancellation of divergences and obtaining a finite 
answer for the matrix element, Aiu, wc perform a very 
strong check on the calculation. Diagram (a) potentially 
has a {aMw)~^ divergence but this vanishes trivially. 
Both (a) and (b) have (aM^y)"^ and {aMw)^"^ diver- 
gences. The (oMvy)""' divergences must cancel between 
the two diagrams, since no such divergences appear in 
any other diagram in Fig. [21 This is confirmed numeri- 
cally, and wc also find that the {aMw)~'^ divergence can- 
cels between these two diagrams. The lattice diagrams 
(c) and (d) contain {aMw)~'^ divergences which cancel 
between them. This is a good check for the derived lattice 
Feynman rules. The derivative diagrams (e)-(g) contain 
modifications of the HW^W~ vertex which restricts the 
Lorentz index on the W bosons at this vertex to be ^ = 0. 
They each clearly contain a {aAlwr)~'^ divergence, and 
we find that the contributions from all three diagrams 
are needed for cancellation. This acts as another strong 
check on the lattice Feynman rules since the contribution 
from (g) is needed to cancel the divergences of (e) and 
(f), thus linking the rule generating diagram (d) with the 
standard lattice rules which give rise to diagrams (a) and 
(b). Indeed, we find numerically that diagrams (e), (f) 
and (g) cancel exactly. 

r dependence — We deduce F{t) by performing a si- 
multaneous fit in (r, aMw) to F{t, aMw) so that the 
continuum extrapolation a — )■ can be done. We ini- 
tially fit the sum of the contributions from diagrams 
(a)-l-(b) and from (c)-|-(d) separately, denoting these 
p(a+b) QT^^ p(c+d) ^ respectively. For each fixed r value 
we compute F{t, aMw) in both cases for values of aMw 
0.003 < aMw < 0.2. We find that F(°+'')(t, aM^y) has 
no discernible aMw dependence in this range at fixed t 
and we find the fit result 

pia+b) ^ 2.088(8)-f0.660(17)r+1.500(3) ^2 - [arcsin 

(8) 

with = 0.25. There is no improvement in fit from in- 
cluding lattice artifact aMw polynomials or logarithms. 
For illustration, we show in Fig. [3] the data for the small- 
est aMw value, aMw = 0.003, and the best fit curve in 

The contribution of the lattice-induced diagrams, (c) 
and (d), does have a dependence on lattice artifacts and 
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FIG. 3. A plot for illustration of F("+'''(r, aMw) for aMw = 
0.003 and the fit from Eq. ^ for F("+'''(r). There are no 
observable lattice artifact effects for any value of aMw used. 
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FIG. 4. A plot for illustration of F^'=+'^^t, aMw) for aMw = 
0.003 and the fit from Eq. Q for F'-''+'^\t). While there is 
a lattice artifact term in the fit, its effect even for the largest 
vale of aMw used is not discernible on the plot being a cor- 
rection of order 1% 



is well-dcscribcd by the function 



as we would expect, and the additional lattice terms in 
give a simple linear behaviour in r which is, 
however, vital to the question in hand namely the value 
of C2. We note also that neither f('^+'')(t) or f('=+'*)(t) 
individually vanish at r = 0. 

A fit to the contributions of all four diagrams (a) - (d) 
gives the final result in the continuum limit of 

F{t) = 1.498(8)+1.000(17)r+l. 500(3) - [arcsinVi^^ ^ 

(10) 

which is in strong agreement with the established di- 
mensional regularization result. We note the fact that 
F{0) = arises directly from our calculation, as it does 
in dimensional regularization, without recourse to the 
Dyson subtraction of other calculations in D = 4 

Conclusions — Physical predictions should not depend 
on the particular way a quantum field theory is regu- 
lated. This paper shows how pcrturbative lattice calcu- 
lations can be used to check the correctness of continuum 
calculations of observables, so long as the questions in- 
volved do not directly involve couplings of gauge fields 
to complex-representation Weyl fermions. Agreement 
between the lattice calculation here of the M^-mediated 
contribution to if —> 77 and the longstanding result ob- 
tained with dimensional regularization should put to rest 
any controversy regarding the Standard Model prediction 
for this decay channel. 

Acknowledgments — We are very grateful to B. C. Al- 
lanach, T. Ishii, and S. Lee for helpful discussions, and 
especially to A. Maharana for collaboration in the early 
stages. This work was supported by STFC under grant 
ST/G000581/1. The calculations for this work were, in 
part, performed on the University of Cambridge HPCs 
as a component of the DiRAC facility jointly funded by 
STFC and the Large Facilities Capital Fund of BIS. 



p(c+d) ^ _o.590(l)+0.340(l)r-0.24(3)(aAfvi/)^ log aMvK , 

(9) 

with = 0.65. The fits indicate that there is no need to 
include a polynomial in aMw and that the lattice artifact 
term turns out to be independent of r. In Fig. |3]we show 
the data for the smallest aMw value, aMw = 0.003, and 
the best fit curve in Eq. For the values of aMw 

used we do not illustrate the lattice artifact dependence 
by plotting the data for other values of aMw since for 
the largest value used the correction is of order 1% and 



unrc solvable on the scale of the graph shown. The well- 
controlled continuum extrapolation is, to a high degree, 
due to the use of the Symanzik improved gauge action. 
The derivative diagrams, (c)-(g), cancel exactly, and give 
no net contribution. For the graphs shown the error bars 
are smaller than the symbol size. 

From Eq. © and ^ we see that F(°+'')(r) contains the 
arcsin term which gives rise to the threshold behaviour, 
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